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Abstract 

A recursion for the joint moments of the external branch lengths for coalescents with multi- 
ple collisions (A-coalescents) is provided. This recursion is used to derive asymptotic expansions 
as the sample size n tends to infinity for the moments of the total external branch length of 
the Bolthausen-Sznitman coalescent. The proof is based on an elementary difference method. 
An alternative differential equation method is developed which can be used to obtain exact 
solutions for the joint moments of the external branch lengths for the Bolthausen-Sznitman 
coalescent. The results for example show that the lengths of two randomly chosen external 
branches are positively correlated for the Bolthausen-Sznitman coalescent, whereas they are 
negatively correlated for the Kingman coalescent provided that n > 4. 
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1 Introduction and main results 

Let II = (IIt)t>o be a coalescent process with multiple collisions (A-coalescent). For fundamental 
information on A-coalescents we refer the reader to [23J and [21]. For n € N :— {1,2,...} we denote 
with lit") = (n^ n) ) t >o the coalescent process restricted to [n] := {1, ...,n}. Note that IB") is 
Markovian with state space £„, the set of all equivalence relations (partitions) on [n]. For £ € £ n we 
write |£| for the number of equivalence classes (blocks) of £. For m 6 {1, . . . , n — 1} let g nm be the 
rate at which the block counting process :— (N^)t>o ■= (|II(™' l |) t >o jumps at its first jump 

time from n to m. It is well known (see, for example, [2UJ Eq. (13)]) that 

n I n-m-l 



771—1 



9nm = A X n - m - l (l-x) m - L A(dx) (1) 



[0,1] 



for all n, m 6 N with m < n. We furthermore introduce the total rates 



n-1 



9n 



m—1 



Ei 1 — (1 — x) n — nx(l — x) 
9nm = / " " 5— ^ " Mdx), T7 6N. (2) 



n-1 



[0,1] 



department de Mathematique, Insitut Galilee, Universite Paris 13, Avenue Jean Baptiste Clement, 93430 Vil- 
letaneuse, France, E-mail: dhersin@math.univ-parisl3.fr 

2 Mathematisches Institut, Eberhard Karls Universitat Tubingen, Auf der Morgenstelle 10, 72076 Tubingen, Ger- 
many, E-mail address: martin.moehle@uni-tuebingen.de 



We are interested in the external branches of the restricted coalescent process IT^™) . More precisely, 
for n G N and i G {1, . . . , n} let 

t im ; := inf{i > : {i} is a singleton block of n|™' ) } (3) 

denote the length of the ith. external branch of the restricted coalescent U\ n > . Our first main result 
(Theorem 1 1.1)) provides a general recursion for the joint moments 

» n (kx, . . . , := E(t£\ • ■ • r„^), j G {1, . . . ,n}, fci, . . . , ^ G No := {0, 1, . . .}, (4) 
of the external branch lengths. The proof of Theorem II .11 is provided in Section [51 

Theorem 1.1 (Recursion for the joint moments of the external branch lengths) 

For all n 6 {2, 3, . . .}, j G {1, . . . , n} and k%, . . . ,kj G N the joint moments fjL n {k\, . . . , kj) := 

^( T n\ ' ' ' T nj) °f th- e lengths T„ r i, . . . ,r n! „ of the external branches of a A-coalescent IL^ satisfy 
the recursion 

l^n{kl, . . . ,kj) 

1 J ™ _1 (m — 1) 

= J ki [l n (k\ , . . . , , fcj 1, , . . . , kj) + \ Prim 7 ^ Mm (^1 ) • • • > kjjt (5) 

where p nm := g nm /9n and g nm and g n are defined via (QJ) and 0j. 
Remarks. 

1. The recursion ([5]) works as follows. Let us call k := ki H + kj the order (or degree) of the 

moment fi n {ki, . . . ,kj). Provided that all the moments of order k — 1 are already computed, 
(|5|) is a recursion on n for the joint moments of order k, which can be solved iteratively. So 
one starts with k = 1 (and hence j = 1), in which case the recursion ([5]) reduces to /i„(l) = 
VPn+El=2Pnm((m-l)/n)/i m (l), n G {2,3, . . .}. Since p 2 (l) = E(r 2)1 ) = l/ 52 = 1/A([0, 1]), 
this recursion determines the moments of order 1 completely. Now choose k = 2 in ([5]) which 
leads to a recursion for the second order moments. Iteratively, one can move to larger values 
of k. 

2. For j = 2 and fei = fea = 1 the recursion © reduces to 

E(r„,ir„, 2 ) = — E(r„.i) + V p nm —, ; j2 E(r m ,iT m , 2 ), n G {2,3,.. .}, (6) 

with initial value E(t2 ! iT2.2) = 2/g|. Provided that E(r nj i) is known, (J6]) is a, recursion on n 
for E(r„ ! iT„ :2 )- 

Note that Theorem 1 1 . 1 1 holds for arbitrary A-coalescents. For particular A-coalescents the recursion 
(J5J) can be used to derive exact solutions and asymptotic expansions for the joint moments of the 
lengths of the external branches. In the following we briefly discuss the star-shaped coalescent 
and the Kingman coalescent. Afterwards we intensively study the external branch lengths of the 
Bolthausen-Sznitman coalescent. For related results on external branches for beta coalescents we 
refer the reader to j7] and [8]. 
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Example. (Star-shaped coalescent) For the star-shaped coalescent, where A is the Dirac mea- 
sure at 1, the time T n of the first jump of IT™' is exponentially distributed with parameter 
g n = 1, n G {2,3,...}. Furthermore, p n \ = 1 and p nm = for n, m G N with 1 < m < n. 
Thus, (O reduces to /7„(fci, . . . ,kj) — Yn=i ki ^n{ki, . . . , fej-i, ki — 1, fcj+i, .. .,kj) with solution 
Unikx, . . . ,kj) = (fci + • • • + fcj)!, which is obviously correct, since r n) j = T„ for all i € {1, . . . , n} 
and, therefore, /!„(&!,...,%) = E(T* 1+ ) = (fci + ••• + %)!, n G {2,3,...}, j G {l,...,n}, 
ki,...,kj G N. 

Example. (Kingman coalescent) For the Kingman coalescent [15], where A is the Dirac measure at 
0, the time T n of the first jump of iT™- 1 is exponentially distributed with parameter g n = n(n — 1)/2, 
7i G {2, 3, . . .}. Furthermore, p n ,n-i = 1, 77, G {2, 3, . . .} and p nm — for all to, n G N with to < n — 1. 
Caliebe et al. [5, Theorem 1] verified that nr n i — > Z in distribution asn-> oo, where Z has density 
a; h-> 8/(2 + a;) 3 , a; > 0. Janson and Kersting [TSJ Theorem 1] showed that the total external branch 
length L e n xternal := £™ =1 T Hji satisfies (l/2)^n/(log7i)(L^ te ™ Q ' - 2) -> jV(0, 1) in distribution as 
n — >• oo. We are instead interested here in the moments of T ni i- The recursion ([3]) for j ' = 1 reduces 
to 

2k n — 2 

A«n(fc) = — —fj, n (k-l) + Mn-iW, n G {2,3,...},fc G N. 

77(77 — lj n 

Rewriting this recursion in terms of a n (k) := 77(77 — l)/i„(fc) yields a„(fc) = 2k fi n (k — 1) +a n _i(fc), 
n G {2, 3, . . .}, fc G N, with solution a n (k) = 2fc^^ =2 /J m (A: — 1). Thus, 

2k " 

M*0 = -p TtZ)M*-1), nG {2,3,...},fcGN. 

n(n — 1 ^-^ 

v ; m=2 

The first two moments are therefore E(r n ,i) = Mn(l) = 2/(n(n — 1)) X)m=2 = V n anc ^ 

Tjji/ 2 \ „ ^ 4 " 2 8(fen - 1) logn 8(7-1) /logn 

77(77 — 1) m n(n — 1) n z n z \ n J 

where 7 » 0.577216 denotes the Euler constant and h n := J^.-, 1/i the 77-th harmonic number, 
n G N. Note that these results are in agreement with those of Caliebe et al. [5j Eq. (2)] and Janson 
and Kersting [18, p. 2205]. For the third moment we obtain 

(3) = 6 y 8 (^» ~ l l = 48 y ' ' 

77(77 — 1) m(m — 1) 77(77 — 1) m(m — 1) 



The last sum simplifies considerably to 



\ h m 1 ^ \ / h m h m 1 

' TO(TO — 1) \ TO — 1 TO TO(TO — 1) 

m-1 v ' m=2 v v ; 

_ hm+1 ^ m ( y ^ 

' TO ^— ' 771 V n 

m— 1 m— 2 v 

n—l-L 1 1 71 1 t , -, 

— /\ -l m+i _ -in n 1 1 

' TO ^— ' TO 77 77 

m— 2 m— 2 
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3 J— * 1 h n 1 

- + > - 1 + — 

2 ^— ' m(m + 1) n n 

m=2 v ; 

3 1 1 h n 1 

- + -1 + - = 1 - — 

2 2 n n n n 



Thus, the third moment of r n ,i is 



. , . . . 48 / fr„\ 48 logn ^ 

n(n — 1) \ n J n z n 6 



For the fourth moment we obtain 



wt 4 n / a\ 8 ,o\ 384 1 - h m /m 

E ( T n.l) = Mn(4) = —, TT >, Mm (3) = -7 TT >. — 7 TT j 

n(n — 1 ' nln— 1) m(m— 1) 

v ; m=2 v ' m=2 v 7 



a formula which does not seem to simplify much further. One may also introduce the generating 
functions fffc(t) := Y^n=2 ^n(k)t n , k £ N, \t\ < 1. For all fc > 2 we have 

oo oo n 

* 2 <?fe(*) - ^^-lK(fc)f = ^2fc^/i m (fc-l)< n 

n— 2 n — 2 m— 2 

oo oo , 

= 2k Mm (A; - l)i m E *"~ m = J3t»fc-i(*). 

m— 2 n—m 

so these generating functions satisfy the recursion 

g k (t) = 2k [ [ 9 *~ 1<yU \ duds, fc>2,0<i<l, 
Vo Jo u 2 (1-m) 



with initial function <?i(i) = ^^2(2 / 'ji)i n = — 2£ — 21og(l — i). Using this recursion, gk(f) can 
be computed iteratively, however, the expressions become quite involved with increasing fc. For 



example, g 2 (t) = St - 4(1 - t)log 2 (l - t) - 8(1 - t)U 2 (t), \t\ < 1, where Li 2 (t) := -J *(log(l - 
x))/xdx = Yl'kLi t k /k 2 denotes the dilogarithm function. In principle higher order moments and as 
well joint moments can be calculated analogously, however the expressions become more and more 
nasty with increasing order. In the following we exemplary derive an exact formula for fi n (l, 1) = 
E(r, i ,ir n 2). The recursion ([5]) for j = 2 and k\ — k 2 — 1 reduces to (see 

. , 2 . , (n-2) 2 , s 8 (n-2)(n-3) 

MM = — Mn 1 + VrM 1 - 1 =^7 TT + 7 Ts V n-l 1,1 , n>2. 

g n \n)2 n*[n — l) n[n — l) 

It is readily checked by induction on n that the solution of this recursion is given by /^(l, 1) = 2 
and 

4(n 2 - 5n + 4/t„) 

Mn(M) = —7 TT!? TyT' rc G {3,4,...}. 

n(n — 1) (n — 2) 

In particular, the asymptotic expansion /x n (l, 1) = 4/n 2 — 4/n 3 + 0((logn)/n 4 ), n — » 00, holds. 
Moreover, Cov(r n ,i,r n ,2) = Mn(l, 1)- (Mn(i)) 2 = 4(n 2 - 5n + 4/i n )/(n(n - l) 2 (n - 2)) - 4/n 2 < for 
all n > 4. Thus, for the Kingman coalescent, the lengths of two randomly chosen external branches 
are (slightly) negatively correlated for all n > 4. We have used the derived formulas to compute the 
following table. 
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n 


M„(l) =E(r„,i) 


fi n (l, 1) = E(T n ,ir ni2 ) 


Cov(r„ i i,r n , 2 ) 


2 
3 
4 
5 
10 
100 


1 

| « 0.666667 

| = 0.4 

I = °- 2 
h = 0.02 


2 

| « 0.444444 
|| k 0.240741 
|| « 0.152222 
Jllg^ 0.038096 
Ri 0.000396 


1 


-jjjjg w -0.009259 
_ V w „ .007778 
-™ « -0-001904 
« -0.000004 


oo 


2 

n 


4 _ 4 + o(i2s_) 


-£ + 0(^) 



Table 1: Covariance of r n ,\ and r n ,2 for the Kingman coalescent 

In the following we focus on the Bolthausen-Sznitman coalescent [4] , where A is the uniform distri- 
bution on [0, 1]. Our second main result (Theorem 11.20 provides asymptotic expansions for all the 
joint moments of the external branch lengths for the Bolthausen-Sznitman coalescent. 

Theorem 1.2 (Expansion for the joint moments of the external branch lengths) 

For the Bolthausen-Sznitman coalescent, the joint moments /i„(fci, . . . , kj) := E(t^\ ■ ■ ■ t^ 3 -) of the 
lengths T n i, . . . , r„ in of the external branches satisfy the asymptotic expansion 



fJ>n(ki, ■ ■ ■ ,kj) 



W.---kj\ 



1 



Kj (k\ , . . . , kj) 
logn 



O 



log 2 n 



j GN,ki,...,kj GNo, (7) 



where the coefficients Kj(ki, . . . , kj), j G N, k\,...,kj G No, are recursively defined via «i(0) := 
0, Kj{ki,...,ki-i,Q,k i+ i,...,kj) := K j - 1 (ki,...,k i -i,k i+ i,...,kj), j G N\ {1}, i G {l,...,j}, 
/«!,..., fci_i, fcj+i, . . . ,kj 6 No, and 



Kj(ki,.. .,kj) := -f J^fc. + - j + l+^«j(fci,.. ■ ,h-iM - l,h + i, . . .,kj)j 

3 ^ i=l i=l ' 

for j,k\,. . . ,kj G N. //ere $ := r'/T denotes the derivative o/logT (a/so called the digamma 
function). 

Remarks. 

1. For j = 1 the recursion for Kj(&i, . . . , fcj) reduces to = fc — 7 + Ki(fc — 1), k G N, 
where 7 = — \&(1) ~ 0.577216 denotes the Euler constant. By induction on k it follows that 
Ki(k) = k(k + l)/2 — kj, k G No, in agreement with [111 Theorem 1.2]. 

2. For ki = ■ ■ ■ = kj = 1 the recursion for Kj{k\, . . . , kj) reduces to «i(l) = 1 + ^(1) + «i(0) = 
1 — 7 and Kj(l, • • • , 1) = ^(j) + 1/.7 + K j-i(l, •••,!) for j G {2, 3, . . .}. An induction on j it 
follows that Kj(l, ■ • • , 1) = (j + l)/ij — j — J7, j G N, where hj := Yli=i V* denotes the jth 
harmonic number, j G N. 



3. For j = 2 and fci = fc 2 = 1 the expansion {?]) has the form 

1 K2(l,l) 



E(r„,ir„ i2 ) = Mn(l,l) 



log 2 n log 3 n 



log n 



00, 



(8) 



with K 2 (l,l) = 5/2 - 27 w 1.345569. In particular, Cov(r„,i, t„ >2 ) = ^ n (l, 1) - (/J„(l)) 2 = 
1/(2 log 3 n) +0(1/ log 4 n), n — > 00. Thus, for the Bolthausen-Sznitman coalescent, r n ,i and 
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i~ n ,2 are asymptotically positively correlated. With some more effort (see Corollary 13.31 and 
the remark thereafter) exact solutions for E(r nj i) and E(r n ir n) 2) are obtained and it follows 
that r n ,i and r n ^ are positively correlated for all n > 2. Note that this positive correlation 
differs substantially from the situation for the Kingman coalescent, where T n> i and T n ^ are 
slightly negatively correlated for all n > 4. 

The following convergence result is a direct consequence of Theorem 11.21 

Corollary 1.3 (Weak limiting behavior of the external branch lengths) 

For the Bolthausen-Sznitman coalescent, (log u)(r n> i, . . . , T n>n , 0,0,...) — > (n, T2, . . .) in distribution 
as n — > oo, where t\,T2,.-. are independent and all exponentially distributed with parameter 1. 

The following result is as well a direct consequence of Theorem ll.2l and concerns the asymptotics of 
the moments of the total external branch length jj^ ternal ~ ^" =1 r ni , of the Bolthausen-Sznitman 
coalescent. 

Corollary 1.4 (Asymptotics of the moments of the total external branch length) 

Fix k £ N. For the Bolthausen-Sznitman coalescent, the kth moment of L<^ ternal has the asymptotic 
expansion 

H(LT ernal ) k ) = 7^( l + r^ + (v^r-))^ ( 9 ) 

log 77 V lo g n \log nJJ 

where Ck '■= (k + l)/ife — k — Ivy. 



Remarks. 

1. Note that c x = 1 - 7, c 2 = 5/2 - 27, c 3 = 13/3 - 37, and c 4 = 77/12 - 47. The expansion © 
quasi coincides with the expansion of the fcth moment of the total branch length L n (see, for 
example, [9j Corollary 4.3]), only the coefficient Ck is of an additive term k smaller than the 
corresponding coefficient mk — (k + l)hk — k-f (see [9]) in the asymptotic expansion of E(L„). 

2. As in the proof of [HI Corollary 4.4] it follows that n~ 1 (log n)U^f ternal — > 1 in probability as 
n 00. The same argument as given in [221 p. 277] yields the asymptotic expansion 

mL e X ternal _ E{L e X ternal ))k) = ^ f + O ( ), k > 2, 

k{k - 1) log L+L n Vlog + nj 

of the centered moments of J J e ^ ternal ; which coincide with those (see [HI Eq. (28)]) of the 
centered moments of the total branch length L n . In particular, 

1 n f 77 \ 

VarfL^"^) = +0 , 77^^. (10) 

^ log n \ log 77 / 



The moments of Jj^ ternal unfortunately do not provide much information on the distributional 
limiting behavior of Jj^ ternal as n — » oo. Nevertheless, the asymptotic expansions of the centered 
moments of L n and Jj^ ternal coincide, which supports (or at least does not contradict) our intuition 
that the distributional limiting behavior of £^ cternal coincides with that of L n . Based on the weak 
convergence result [9] Theorem 5.2] for L n we therefore state the same convergence result for 

^external &g & con j e cture. 
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Conjecture 1.5 (Weak convergence of the total external branch length) 

For the Bolthausen-Sznitman coalescent, 

1 2 

^B_» L «rter«oI _ log n - log log 71 -¥ L (11) 

n 

in distribution as n — > oo, where L is a 1-stable random variable with characteristic function t t— > 
exp(it\og\t\ -7r|i|/2), t 6 R. 

Remark. The same scaling and limiting behavior as in (fTTj) is known for the number of cuts to 
isolate the root of a random recursive tree ([TU], [IB])- Essentially the same scaling and convergence 
result has been obtained for random records and cuttings in binary search trees by Holmgren [IH 
Theorem 1.1] and more generally in split trees (Holmgren [T2] Theorem 1.1] and [T|3 Theorem 1.1]) 
introduced by Devroye [BJ . The logarithmic height of the involved trees seems to be one of the main 
sources for the occurrence of such scalings and of the 1-stable limiting law. 



2 Proofs 

Proof, (of Theorem 11.11) Let T n denote the time of the first jump of the block counting process 
N^ n > and let I n denote the state of the block counting process after its first jump. Note that 
T n is exponentially distributed with parameter g n and that P(/ n = m ) = Pnm '■= 9nm/gn for all 
to G {1, . . . , n — 1}. Moreover, T n and I n are independent. We verify ((SJ by induction on the number 
j of involved external branches. 

Step 1. The case j = 1: 
We have to verify that 

k n ~ 1 — 1 
E(r^) = _E(t*i 1 )+ y> nm - E(t* iX ), 7i e {2,3,...},* GN. (12) 

Note that (|12[) was already verified in [11J . however we provide a proof which turns out to be 
generalizable to an arbitrary number j of external branches. We verify (1121) by induction on fc G N. 
Conditional on I n = m, the external branch 1 is involved in the first collision event with probability 
qi := (n — m + l)/n and it is not involved in the first collision event with complementary probability 
qo := (to — l)/n. Thus, 

E« 1 |7 n =m) = qi E{T^) + q E{{T n + T mA ) k ), 



where r m .i is independent of T n , Binomial expansion yields 

E(r n fe x | J„ = to) = qi E(T 7 k ;) + qo J2 (f)E(2$)E(7£?) = E(T*) + <7o £ ( " )E(^)E(r^). 



i=0 



Multiplication with p nm = V(I n = to) and summation over all to G {1, . . . , n — 1} leads to 
E«0 = E(T^ + J2(%(K)J2Pnrn^^E(r^) 



n 

i=0 v 7 m=2 



fc— 1 \ n— 1 n — 1 1 

E^) + ^ EK) ^ P „ m — E(i ! ) + ^ Pnm — E(r^). (13) 

j=l v 7 m=2 m=2 
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For k = 1 this equation reduces to 

1 n-1 -. 

1 » to — 1 . 
E(t„i) = h > »„ m E(r mi i). 



<7n Jl 

y ™ m=2 



Thus, (fT2"j) holds for fc = 1. The induction step from {1, . . . , k — 1} to fc (> 2) works as follows. By 
induction it is allowed to replace the first sum over m in (fT3|) by E(t„^ j ) ~ {k — 1 E(r^ t_1 ). 
We therefore obtain 



E(r* tl ) - E(T*) + W^Wl)(E(^ 

i=l ^ ' ^™ m=2 71 

- E (^E(^)E(r„V) - E (^W^E^ 1 ) + g p^Z^E^) 

* 1 V / A 1 V / " n <m — 9 



z— 1 x 7 m— 2 



= feE(T„)E(r^ 1 ) + ( ■) ECW<?) - £ ( k V^^-i^E^*) 



2=2 

n— 1 ^ 

m=2 

= AE^ + l^n™— E^), 

since, for all i G {2, . . . , fc}, 

' fc \^^-n fc -* + 1 _ AV ! f k \ {i-l)\k-i + l 



vV V*- 1 / 5n Wffn V 1 " 1 / ffn ffn 

This finishes the induction and completes the proof of (fT2|) . Thus, ([5|) holds for j = 1. 



The induction step from to {1, . . . , j — 1} to j (> 2) for general number j of involved external 
branches is somewhat technical. We therefore consider next the case j = 2. 

Step 2. The case j = 2. 

The proof works again by induction on the order fc := k\ + ki- We thus first have to verify (J6j) , 
which is the particular case j — 2 and fci = &2 = 1 of the general recursion (|5|) . 
Conditional on I n = to, the two external branches 1 and 2 are both involved in the first collision event 
with probability qn :— (n — m+ 1)2/(^)21 the external branch 1 is involved and the external branch 
2 is not involved in the first collision event with probability q±o := (n — m+l)(m — l)/(n)2, and both 
external branches are not involved in the first collision event with probability qoo := (to — 1)2/(^)2- 
Note that qiQ — goij that gn + gio + qoi + goo = 1 arid that q±o + goo = [m — l)/n. Thus, 

E(r ni ir n ,2 I I n = m) = gnE(T^) + q w E(T n (T n + T m . 2 )) 

+goiE((T„ + T m ,i)T n ) + gooE((T„ + r TO ,i)(T„ + T m , 2 )) 

= E(T*) + 2E(T n )E(r m , 1 )(g 10 + goo) + 9ooE( 

2 2„, ,m-l . Jm-l) 2 

= — + — E(r m 1) + E r m ,ir m ,2 V , v ■ 

3n 5n n (n)2 
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Multiplication of both sides with p nrn = P(/„ = m) and summation over all m £ {1, . . . , n — 1} 
yields 

2/1 ^-4 to — 1 \ ^-4 (m — 1)2 . 

E(T n> lT n ,2j = h > Pnrn E(r m ij + > p nm -—■ E(T TO ,iT m , 2 ) 

9n\9n ^ n J ^ (n) a 

which is ([6]). Concerning the induction step from {1, . . . , ki + fc 2 — 1} to k := k\ + fc 2 G {3, 4, . . .} 
we have to verify that 

¥ ( T ki k 2 n _ ^l Ff fex-1 fc 2 n , fex fe2-l^ , V n ( TO Z 1 ) 2 Fr T fcl -r* 2 ^ 

9n 5« ^ (n) 2 

We have 

E ( T n'i T n% I 4 - m) = 9ll E(T^+ fe2 ) + g iE((T„ + r m4 ) fel T^) 

+ (Zl0 E(T„ fel (T„ + r m , 2 ) fe2 ) + g oE((T„ + r mA ) kl (T„ + r m , 2 ) fe2 ). 

Binomial expansion together with the fact that T n is independent of (T mi i,T mj2 ) leads to 

fcx 



E«\r5|7„ = m) = gil E(T*+«*) + q m £ (^W^E^ 1 ) 

£ Me(T^)E(t^) 



i 2 =0 
fex k 2 

-qoo 

ix=0 i 2 =0 

fex-1 



E(r* 1+fe *) + (goi + <Zoo) £ (f)E(TWfe)E(^r 



ix=0 

fe 2 -l 



'1 



(9io + goo) £ (^W^E^ 2 ) 
» 2 =o ^ 2 ^ 

EE f? 1 ) ^ 2 V(^ + -)E(^r-^ 2 --). 



-9oo 

tx=0 i 2 =0 



«1/ V«2 



Note that qio = qoi and that goi +P00 = (™ ~ l)/ n - Multiplication with p„ m = P(/„ = to) and 
summation over all m £ {1, . . . , n — 1} leads to 

fex — 1 /, v. n — 1 

E«\t5) = E(T^ +fc2 )+E( 1 EK^)^^— E(r^) 

»l=Q ^ Y ' m=2 
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k-2 — 1 /7 \ 71—1 



i 2 =0 ^ 2 ' m=2 '' 



ii — *2—0 rn—2 

By (HU), the first two sums over m are equal to E(r^f ll ) - (A* -i^^E^f <I_1 ) and E(r^ 2 _42 ) - 
(fe — i2)5^ 1 E(r^ 2 2 ~ l2_1 ) respectively. Moreover, by induction it is allowed to replace the third (last) 
sum over m (except for the case h = i 2 = 0) by E(r£ l f ^r^ 2 ) - (fci - H^E^f il_1 r^ 12 ) - 
(fe - i2)g r ; 1 E(v 1 r il r^ 2 2 - 12 " 1 ). We therefore obtain 

E«W fe2 2 ) = E(T^+ fe2 )+E f fcl )E(T^+ fc2 )(E(r r ^ 1 )-^^IE(r5- 11 - 1 )) 



ii=0 

fe 2 -l 



12=0 



E ( fc2 )E(T^+-)(E(r„V 2 ) - -i-^E^r 2 " 1 )) 



E (£)(*W-)x 



'm?:: 1 , 'v::-, ■-) mrw "r:-.-, - j — - iir^ 1 , 1i t*\ 12 1 ' 



i V d 1/2 E(r fcl r fe2 1 

m=2 ^ > 2 

The last sum is already what we need. For the other sums the main task is to reorder all the 
summands in terms of the moments of r„ i and r n ^ 2 . The coefficient in front of E(T^ 1 1 ~ 1 T r j 2 2 ) is 

( fe Q 2 )E(T^+°) = ki/g„ and the coefficient in front of E^T^f 1 ) is ft 1 ) ( fe ] 2 )E(T°+ 1 ) = k 2 /g n as 
required. It remains to verify that all the other coefficients are equal to zero. Careful investigation 
shows that the coefficient in front of E(r T j 1 1 ~ n ) (with %\ < kx) is 



E(T£ 1+fe2 ) 



( h )E(T^-^) kl ~ {H - 1) -( kl )( k2 W^* 2 " 1 )) 
\ti-lj " g n \iij\k 2 -lj 

ki\(h + k 2 )\ ( h \ (ii + fc 2 - 1)! k% - h + 1 fkA (ii + fe 2 -l)!l 



+ (fc2 -l h fc2 - (fa - 1) 



'1/ g n \tl J-/ <?„ <M 5n i/n 



Al\ (ti + fa-1)! ^, , v 



9)i 

In the same manner it follows that the coefficient in front of E(t^ 2 2 -12 ) (with i 2 < k 2 ) vanishes. The 



coefficient in front of E^Y"* 1 ^ - * 2 ) (with 1 < h < kx and 1 < i 2 < k 2 ) is 



'fcA fk 2 
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ki\ (k 2 \ (h + 



1) 



HJ\i2j gl 1 ^ 2 



ki \ («i + i% — 1)! k\ +ti + 1 (k\\( k 2 \ (i% + i 2 — 1)! — *2 + 1 



The last coefficient to consider is that in front of E(r° 2 ) (= 1), which is 

E(2*+*»)-( ^ )E(T r ^-^) kl ~ {kl - 1] - ( k2 ^^MljMM 

\fcl - 1/ 3" \ fc 2 ~ 1/ 9n 

(fci + fc 2 )! (fa + k 2 - 1)! 1 (fci + fc 2 - 1)! 1 

fc 1+ fc 2 fc l fc 1+ fc 2 _l fc 2 fc +fc2 _ a U. 

9n gn gn <M 

This finishes the induction and the proof of for j = 2 is complete. Note that we have used (fT"2"|) . 
so in order to verify ([3]) for j = 2, we have used ([3]) for j = 1. 

Step 3. The general case j > 2: 

The proof of the induction step for general number j (> 2) of involved external branches works 
essentially the same (induction over the order k := kx + • • • + kj of the moments), however, the 
details become more involved and a bit nasty to write down. We leave the details to the interested 
reader. □ 

Proof, (of Theorem II. 2 1) For the Bolthausen-Sznitman coalescent, g nm = n/((n — m)(n — m+ 1)), 
to e {1, . . . , n — 1} and g n = n — 1, n € N. Thus, the random state of the block counting process 
after its first jump has distribution 

77 

Pnm = P(/„ = to) = — — , n 6N,m G {1, ...,n- 1}. 

(n — l)(n — m)(n — m + 1) 

We verify by induction on the sum fc := fci + • • • + fcj € No- Obviously holds for k = 0. In 
order to verify the induction step from k — 1 to k € N fix fci , . . . , kj £ No with fci + • • • + kj = k and 
for n 6 N define a n := /_i„(fci, . . . , kj) for convenience. We apply the difference method used in [17] 
to the recursion {5j , which is of the form 

E(to — l)j 
Pnm ~, 7 
\ n )j 



m— 1 

n— 1 / -.n 

71 (TO — 1) 



E 



' (n — l)(n — m)(n — to + 1) (n)j 



J Or. 



E 



(m - l)j 



(n - l)(n)j ^ (77 - to) (77 - to + 1) 
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where 

1 j 

Qn := ^ k{fl n (fcl i ■ • ■ j ki — 1 , k{ 1 , fcj-f 1 , . . . , kj / 



1_ -yy fcj! ■ ■ ■ kjl / Kj(fci,...,fci_i,fcj ^ l,fci + i,...,fcj) ^ 

1 fcl! • • • kjl y—^ ( Kj(kl, . . . , fcj — 1, fcf+1, . . . , fcj) 



' ' ' ^j 1 / ^ . ^jl^ll ■ ■ ■ , ft-j — i-, ■ ■ ■ , l^j ) Q 



g n log n ^ V lo S"- \log n 

From l/g„ = l/(n - 1) = ^(l/n)' = 1/n + 0(l/n 2 ) it follows that 

_ jfci! • • -fcj! | fci!---A;j!X)i = iKj(fei)---)Ai-i)Ai - ••-,%) , ^/ 1 



nlog n nlog n \nlog + n 

The difference b n :— a n — k\ \ ■ ■ ■ kjl / log* n satisfies the recursion 
hi- ■ -k^l 



log n 



(m-l)j ( fci!---fcj!\ ki\---kj\ 



(n- l){n)j ^ (n-m)(n-m+l) \ log* m 7 log*n 

7i — m -\- 



(n - l)(n)j (« - m)(ri - m + 1) 



with 

4 



^— J (m— l)j fci!---fcj! 



(n - l)(n)j ^ 2 (n - m)(n - m + 1) log* to log* n 
Corollary 14.21 in the appendix yields 



n-1 



ri jn?- 1 /, , T ... + ^ ^ 



_ 2 (n — m)(n — to + 1) log* to log* n log* 1 n \ " 2 /log*n Vlog* +1 n 

and hence 



n (to — l)j 



(n - l)(n)j ^ 2 (n - m)(n - to + 1) log* 



n 3 ni +1 \ n> 
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X 



- Jprr + (* + 5*0) ~ j0 + 1)/2)-^ + o(-n^) 

\ log rt log n log n v log n ' / 



-4r^+(fc + i*(i)-i + i)-V + ( 1 l+i 

log n n log n \ / n log n \ n log n 
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We therefore obtain 



fci!-- - kj\( y^«j(fei, . . .,ki-x,ki - l,k i+ i, ...,kj) + k + j9(j) - j + 1 j — — !- 



i=i 

(>[ -An-) = 3K—^- + o( ' 



n log fc n 



v nlog fc+1 n/ nlog fc n \wlog fc+1 n / 

with 

Jf := fc x ! • • -kj\- (k + j^(j) - j + 1 + ^2 Kjih,. . . , - 1, fej+i, . . . J 

^ ^ i=l ' 

= k%\ ■ ■ ■ kj\Kj(k%, . . . , kj ) 

by the recursive definition of the coefficients Kj(k±, . . . ,kj). We now proceed with the recursion for 
b n in the same manner. The difference c n :— b n — K j log t+1 n satisfies the recursion 

c - b - K 



log fe+1 



n 

n-l 



ii (m — 1), / K \ K 

= gn + 7I Tw3 E ' 



m=2 
n-l 



(n-l)(n)j^' 2 (n-m)(n-m+l) V ™ log fc+1 m/ \og k+1 n 
m. — ? 

m=2 



9n + (^T)H~ f-, (n - m)(n - m + 1) 



with 



n \ (m — l)j K 



(n - l)(n)j ^ ( n - m)(n - ro + 1) log fc+1 m log fc+1 
Corollary 14.21 yields 

n-l 



m=2 

such that we obtain 



y ( TO - j_j = n " j ft " , _ 

Z- "' (n — m)(n — m + 1) log fe+1 m log fc+1 n log fe n \log fc+1 n 



C = (jK-M-^ + oi—^r-) = Of 1 



ifc ■ - \ 1 fc+1 / I 1 fe+1 

n log n \n log nj \n log n 

By Lemma 14.31 applied with k replaced by fc + 1, it follows that c n = 0(1/ log fc+2 n), and the 
induction step is established. □ 

Proof, (of Corollary II .3j) Theorem 11.21 clearly implies that, for j £ N and fc x ,...,fcj G No, 
E((r, M logn)^.--(r nj logn) fc O = (logn) fcl +-+**/i n (fei, . . . , kj) -> /,,!••• /,;,,! = Ifa* 1 • • • r*' ) as 
n ->• co. For alii e {!,... ,j\ and allO < 6 < 1 we have E^o(^7 r! ) E ( r D = E^Lo 6|r = 1 /( 1 ~ d ) < 
oo. Therefore (see [2], Theorems 30.1 and 30.2 for the one-dimensional case and Problem 30.6 on 
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p. 398 for the multi-dimensional case) the above convergence of moments implies the convergence 
(log n)(r n) i, . . . , T n j) — > (n, . . . , tj) in distribution asn^oo for each j 6 N. The convergence of all 
these j-dimensional distributions is already equivalent (see Billingsley [31 p. 19]) to the convergence 
of the full processes (log ra)(r n) i, . . . , T ni „, 0, 0, . . .) — ► (ri, T2, . . .) in distribution as n — > oo. □ 

Proof, (of Corollarv ll.4p The total external branch length jj^ ternal satisfies (see, for example, [2TJ 
p. 2165] 

mL e,ternal )k) = £ M £ ^ Mfa, . . . , fej), U S {2, 3, . . .}, fe S N. 

j— 1 V^/ fc 1 ,...,fe- eN ^* J* 
fc j • ■ - -j-fcj =fc 

By Theorem rOI it follows that 

From this representation it follows immediately that asymptotically all the summands with in- 
dices j < fc can be disregarded, so only the summand with index j = k is asymptotically of any 
importance. Thus, we obtain 

mL extemal )k) = (A Jd / + Kk (l, . . . ,1) + Q { 1 



k) log fc n \ log n \ log 2 n 

n k ( «*(!,..., 1) 



1+ »v » . / + Q 



log k n\ ' logn ' Vlog 2 n, 
where «;&(!, . . . , 1) = (fc + l)/ife — fc — fc7 by the second remark after Theorem 1 1.2 1 □ 



3 Differential equations approach 

The expansions of Theorem 11.21 in Section [1] are essentially based on the difference method used in 
the proof of Theorem 11.21 An alternative method based on generating functions is now provided. 
This approach yields for example exact expressions for E(T ni iT n ,2) m terms of Stirling numbers. For 
j e N and fc = (fci, . . . , kj) £ N J define the generating function 

oo 
n=j 

Note that f[ l) (0) = for I e {1, . . . ,j - 2} and that /^ _1) (0) = (j - ■ • • rfy). We shall 

see soon, that the functions fk are analytic in the domain D := {s <E C : \s\ < 1}. The functions 
fk even have analytic continuations on C \ [l,oo), but we will not use these continuations in our 
approach. We start with the following lemma, which is well known (see [111 Lemma 3.1, Eq. (3.3)]). 

Lemma 3.1 For all z € D, fi(z) = / j- — - — 7- y-dt. 
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Proof. Define the auxiliary function 

.W:=f;^" = i-Mi-,) + ^lz£), ,eD. 

71=1 V ' 

For zeflwe have 

oo 

zf[(z) = ^(n-llEK,)^ 1 



=2 

oo 



n=2 

OO OO 71 — 1 1 

oo n— 1 i 



OO / 

Z 



1 . 2 .Ei ra E(«-'-w^"-'- 1 

(=1 V ; n=Z+l 



oo 

Z 



EttttttE^- 1 )^)^- 1 



1 - z f-' i(i + 1) ^ 

Z=l v ' m=l 

= ^ +a (z)z/{(z). 
1 — z 

Thus z(l — a(z))f[(z) = z/(l — z) or, cquivalently, 

1 z 
f[{z) = (l-z){l-a{z)) = (l-z)2(-log(l-z))' ZGA 

since 1 — a{z) = (1 — z)(— log(l — z))/z. The result follows by integration, since /i(0) = 0. □ 

In the following a solution for /(^i) is presented. As a corollary, an exact formula for E(t„ i iT„ 2) is 
provided. 

Lemma 3.2 The function f := /(i,i) satisfies the differential equation 

(l-z)(-log(l-z))/"(z) + /'(z) = (1 _ z)2( _^ og(1 _ z)) , zeD. (14) 
with initial conditions /(0) = and /'(0) = 2. T/ie solution f satisfies 

™ = ^sr—) [ (i-w-Mi-o) * z e D \ (0} - (15) 
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Proof. For n G {2, 3, . . .} we write a n :— E(T ni iTn, 2) for convenience. Note that f(z) = 
S^ =2 a„z™- 1 . In particular, /(0) = and f (0) = a 2 = E(t 24 t 2 , 2 ) = E(T 2 2 ) = 2/.g| = 2. For 
z € D we have 



z 2 /"(z) + z/'(z) = £(n-l) 2 a„* r ' 



n=2 

00 _ n— 1 



„n-l 



£>-l) 8 (— E(t„,i)+ E- , 

K 9n ^ (n) 2 

2$> - lJE^i)^- 1 + £ E ( (W ,'. 1)(W ~ 2) ^ ^^ 1 



n=2 



^— ' ^— ' (n — m)(n — m + 1) 

n=2 m=l v ;v y 



OO 1 OO 



f^Z(Z + l) z 

i=l v ; n= 



i+1 

1 00 



= 2z/;(z)+^— -^(m-l)(m-2) 

;=i ' m=l 

= 2z/ 1 (z) + a(z)z 2 /"(z), 



or, equivalently, 



z 2 (l-a(z))/"(z)+z/'(z) = 2zm = (1 _ z)2( 2 f Qg(1 _, )r 

which yields the differential equation (fH)) , since 1 — a(z) = (1 — z){— log(l — z))j z. Division of (fTfj) 
by 1 — z shows that 

log(l - z)/"(z) + -^/'(z) = -£(-log(l _*)/'(*)) 



(1 - z) 3 (-log(l - z)) bv ' J y ' l-z J K ' dz 

and integration together with the initial condition /'(0) = 2 yields (|T5|) . □ 

Corollary 3.3 (Exact formula for E(r„ iT n 2 )) 

-Fia; n G {2, 3, . . .}. -For i/ie Bolthausen-Sznitman coalescent, 

2 2 fe - l 

Efo.xrn.a) = 7— TvE^^ 3 ^- 2 ^- 1 )' ( 16 ) 

v fe=i 

where the s(n, fc) denote the absolute Stirling numbers of the first kind. 

Remark. Together with the exact formula E(r n ,i) = ((n .— 1)!) _1 Y^k=i s ( n ~ 1) k)/k for the mean of 
T n ,i (see, for example, Proposition 1.2 of [11]) it can be checked that Cov(r n ,i, r„, 2 ) = E(t„,it„ i2 ) — 
(E(t„,i)) 2 > for all n > 2. Thus, for all n > 2, t„ i and r„ i2 are positively correlated. We have 
used the exact formulas for E(t Mi i) and E(r„,iT n 2 ) to compute the entries of the following table. 
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n I E(r n ,i) 



E(r n ,ir n>2 ) 



Cov(r n> i,r n ,2) 



2 
3 

4 
5 

10 
100 
200 
300 
400 
500 
1000 



1 

| = 0.75 

H ps 0.638889 

§ fa 0.572917 

yb 

ps 0.431647 
ps 0.228368 
ps 0.198537 
ps 0.184283 
ps 0.175300 
ps 0.168891 
ps 0.151582 



| = 0.75 
2| ps 0.509259 
ffl ps 0.397569 
ps 0.215119 
ps 0.057067 
ps 0.042758 
ps 0.036676 
ps 0.033092 
ps 0.030652 
ps 0.024546 



^ = 0.1875 
^ ps 0.101080 

1296 

^ ps 0.069336 

0.028800 
ps 0.004915 
ps 0.003341 
ps 0.002716 
ps 0.002362 
ps 0.002128 
ps 0.001568 



I 5/2-2 7 n , 1 x 
log 3 n + ^Ee*^ 



logn log 2 n ^ log' 5 n > 



1 

log-'' n 



Ml 



2 log a 71 1 ^ V log 4 



Table 2: Covariance of r n! i and r„,2 for the Bolthausen-Sznitman coalescent 

Proof, (of Corollary I3.3|) The substitution u = — log(l — t) in the integral on the right hand side 
in (|15|) yields 



/'(*) 



-log(l- 


*) 


2 




-log(l- 


z) 


2 




-log(l- 


z) 


2 





- log(l-z) e 2u _ p v 



du 



■log(l-*) 1 r ~ (2u jfc ~ fc^ 

v fc=0 fe=0 
°° 9 fe - 1 /"-'ogf 1 - 2 ) 



log( i_ z) ^ 



^ 2 fc -l(-log(l-z)) fc 



fe=i 



2* - 1 
kkl 



(-log(l-z)) 



fe-i 



This series expansion of /' in particular shows that / has an analytic continuation on C \ [1, oo). 
From (see [TJ p. 824]) (-log(l - z)) k /k\ = £,°l k ^/i^s(i,k) we conclude that 



/'(*) = 2^ 



2 fe - 1 



__ z i i±i 2 fe - 1 



E ?r-(^-i) - ^E 



s(i, fc — 1). 



fc=l i=fc-l i=0 fc=l 

Thus, using the notation [z l ]f'(z) for the coefficient in front of z 1 in the series expansion of f'(z), 
(i + l)a i+2 = [zV'(z) = 



v. * — ' k 2 

k=l 



or, equivalently, 



2 l±i 2 fc - 1 



Oi+2 



+ 



nyE^^^fc- 1 )- 
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It remains to substitute n = i + 2. □ 
Final remarks and open problems. 

1. Based on the solution (|T5|) for /n n it is possible to derive the asymptotic expansion © of 
E(T ?li iT nj 2) alternatively via expansions of the functions and integrals occurring in (|15p . We 
leave the details to the reader. 

2. Without to much effort it should be possible to derive differential equations in analogy to (TT4"1) 
for fk for arbitrary k — (fci, . .. ,kj) and to gain information on (the coefficients of the series 
expansion of) fk from these differential equations (exact solutions and asymptotic expansions 
of the coefficients via singularity analysis). 

4 Appendix 

Basic expansions (Lemma 14.11 and Corollary 14.21) are provided which are needed in the proof of 
Theorem II .21 We furthermore provide a result (see Lemma l4.3p . which may serve as a stopping rule 
for difference procedures of the form used in the proof of Theorem 11.21 

Lemma 4.1 For all k 6 N, as n — > oo, 

n-2 

'V* n o( n ~\ (17) 

j~[ \og k {n-l) ~ \og k n + Mog^J' 

n ~ 2 11 1 ,1 

1 1 _ 1 , 7 , „ 1 



I 

n-2 



^ I log fc (n - I) log*' 1 n log k n \\og k+1 J 



l + Hog k (n-l) log* _1 n + log k n + °Mog fc+1 n )' ^ 



" ~ 1 1 1 /•• 1 ../ 1 



and 

^ l(l + l)\og k {n-l) ~ log k n + nlog k n + ^nlog k+x n)' ^ 

Proof. Eq. p7|) is a special case of Panholzer's summation formula (see [221 Lemma 4.1, Eq. (16)]) 
and probably as well a consequence of Euler's summation formula (see, for example, |12[ p. 469]). 
We prefer here to verify (fTTj) via a direct decomposition method, which will turn out to work as 
well for (|T5J), and (j2"0")) . We basically split sums of the form ^"=1 • • • ^ w0 P ar ts ■ ■ • 

and X)«=a +i ■ • •' an< ^ nan dle these two parts separately. We will work with the sequence (a„)„ e N 
defined via a± :— 1 and a n := n — |/i/log fe+1 n \ f° r n>2. We will furthermore need the constants 
Cfc := 1/ log fe 2 > 0, k 6 N. In order to verify (fT7f it is sufficient to show that 



n-2 1 

y ( I L_) = f " ) 

^Vi og fc ( n _/) log fc n / \\og k+1 n<> 
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The sequence (a n ) n ^n is chosen such that 



o < E —) < E -r 1 — 

^ W(" - I) log* J ~ l= ^ +1 \og\n I) 

n-2 

< E Ck - Ck(n-a„) = 0[- — 7-n— 

— MOtr 77 

2=a„ + l 1U & ' l 



Moreover, 



^ log fe 77-log fc (77-0 



§(iog fe (77-0 l0g fc n) ^ log fc 77log fc (77-0 



< e 7 V(log fc n-log fe (77-n) 

- log fc nlog fc (77-a„) £f 

1 n— 1 
log 77 l=1 

since log(n — a n ) ~ logn as 7i — > oo. Thus, in order to verify (|17p. it remains to show that 

n-1 

J] (log* n -log' (77-0) = 0(nlog* -1 n). (21) 
i=i 

Applying the formula b k - a k = b k - (a - b + b) k = - J2j =1 ( k ) (a - b) j b k - j with a := log(77 - Z) and 
b := log 77 yields 

n-1 n-i fe /A 

^(log fe 77-log fc (T7-0) = -EE log fc - J r7log-''(l-V«) 



1=1 Z=l j"=l 



fe ^jA n— 1 



= E :)log fe - J 'r 1 ^(-log'(l-V«)). 

3=1 (=1 

Since for each j 6 {1, . . . , k}, 

, n-i ~i 

~yv_ log* (i _//„)) -> / (-lo^ci-x))^ = (-iy +1 j\, 
71 i=i J ° 

it follows that (f2"Tj) holds, and (fTT)) is established. 

Let us now verify (fT8|) . We essentially apply the same decomposition method. Note that ^^1=1 V 
logn + 7 + 0(l/n). Thus, (fT8|) holds if we can verify that 

yl( I L-) = of^_ 

^ l\\og k (n-l) log* 77' Mog fc+1 n 
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The sequence (a ra )„ e N is chosen such that 

^ E y d */ n ~r^) - Cfe E T = °(r^^)- 

Moreover, 



llog fc ?i-log fc (n-Z) 
^ i Vlog fc (n - I) \og k nJ log fe nlog fe (n - V) 



yl( I L_) = v 

^ I \\os k (n-l) loff fc n/ fr* 



< 



1 ^ log fc n-log fc (n-Z) 



E 



log nlog (n - a„) J=1 i 
1 ^ log fc n-log fc (n-Q 

1 2fc 7 ' 



since log(n — a n ) ~ logn as n — > cxd. Thus, it remains to verify that 

' log fe n-log fc (n-0 0/1 fc _ x 



^ 1US = 0(log fe -V). (22) 



Applying the formula b k - a k = b k - (a - b + b) k = - £^ =1 Q)(a - b) j b k -J to a := logfn - I) and 
b := log n yields 

Elog n — login — I) \^lv^/' ! l, fc-i , , ; x 

— 7^ = -EyE (i l0 S k J nhg J (l-l/n) 

1=1 1=1 j=i ^ J ' 



3=1 (=1 



Since for each j 6 {1, . . . , fc}, 

1=1 ' ^ 21 

it follows that (|22p holds and (fl8|) is established. Eq. (|19l) is shown similarly. Let us finally come to 
([2"0jl. Note first that 

^ 2 -log(l-//n) l^ -bg(l-l/n) _ f 1 - 1 /" -log(f_x) 

Kl + l) n ^ (i/n)(i + l)/n " i 1/n r« + 

(f - .T)log(f-a;)i 1 - 1 /« 



L X Jl/n 



and, therefore 
n— 2 



y—(— 



-log(l - Z/n) 



-{ 1(1 + 1) Vlog fe n log fe+1 n 
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1 ^ /l 1 \ k \^ ~ log(l — l/n) 



1 \ k ^ -log(l- l/n) 



log fc nV n-l) + log fe+1 n^ i(/ + 1) 

V- f 1 - - + of \)) + -4r- (— + °C-) 



log n v n Vn 2 // log + n 
1 , 1 



log K n n\og K n \n\og k+1 n 



Thus, in order to verify (|20)) , it suffices to verify that 

y^ 2 1 / 1 1 h -h 

1(1 + 1) \\og k (n-l) \og k n log k+1 n J "V n l og fc+1 n 

The function f n k : [0,n) — > R, defined via 

ftfxl ._ 1 1 -logfl-s/n) 

log (n — a; J log n log n 

has derivative ^ ^ 

= ^(log fc+1 (n-x) _ io?^) " ° 
and satisfies f n k{0) = 0. Thus, f n k > and it follows that (a n ) ne i!j is chosen such that 

„ < g 'f_J___L_ _,r«i-'/») 



< 



^ +i Z(/ + l)Vl O g fe (n-0 log fe n log fc+1 n 
V 1 1 < (I 1 



*4r + i^ +1 ) lo s>-o 

/ 1 1 \ n — a rl — 2 

Cfc — = Cfc 



a„ + 1 n-l/ (a„ + l)(n — 1) 



n/log fc+1 n / 1 



c fc ^^ = O 



l fc+i 
n log n 



Moreover, 



£iaTi)(i 



fe -log(l-Z/n) 



^ 1(1 + 1) Vlog fe (n - I) log fc n log^ 1 n 

1 log fe n - \og k (n - I) - fe(-log(l - l/n))\og k (n - l)/\ogn 
,^ KlTT) log fe nlog fc (n-/) 

< 1 y> log fe n - \og k {n - I) + felog(l - l/n) log fc (n - I)/ logn 



log fe nlog fc (n- a „)tr + 1) 
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— Y 



log te n - log fc (n - I) + fclog(l - l/n) \og k (n - l)/\ogn 



\og 2k nf^ 1(1 + 1) 

since log(n — a n ) ~ logn as n — > oo. Thus, it remains to verify that 

^ log fc n - \og k (n ~ I) + fclog(l - l/n) \og k (n - I)/ \ogn _ n /log^n\ . . 

£j Z(Z + 1) ~ V n J' 

Applying the formula b k - a k + k(a - b)a k /b = (1 - k(a - b)/b)(b k - a k ) + k(a - b)^- 1 = (l-k(a- 

&)/&)(- E?=i Qia-byb^+Ha-b)^- 1 = E-i 1 (fcO-i)-(i))(°- & ) i6 *~ j witha := lo s("-0 

and b :— \ogn yields 

< E ^ " ~~ ~~ ^ + k ^ og< ^ 1 ~ l / n )\°g k ( n ~ 0/ lo S n 



i„„fc 



< 



E 



\og k n- \og k (n - I) + fclog(l - Z/n) log fc (" - l)/\ogn 



1(1 + 1) 



1 = 1 

n-1 ^ fe+1 

i=i ioTT)§v"Vi-i; 

Since for each j € {2, . . . , k + 1}, 

^ log*(l-i/n) _ l^ 1 log^l-i/n) /- 1 Wq-a;) 

^ 1(1 + 1) ~ n^(l/n)((/ + l)/n) ~> 7 * 2 ' 

it follows that the expression on the left hand side in (|23l) is even equal to 0((log fc ~ 2 n)/n). In 
particular, (|2"5|) holds, and, hence, (|2U)) is established. □ 

Corollary 4.2 For aZZ j, k 6 N, as n ^ oo, 

(to - l)j 



E 



(n — to) (n — m + 1) log to 

_ i J J— ;— + (fe + J*C?) — -) r h 0( EXi— )• 



log n log n ^ 2 / fog n ^log n 

Proof. We have 



^ (n - m)(n - to + 1) log to ^ Z(Z + 1) log (n - Z) 
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Plugging in 

(n-l-l)j = JJ(n-(I + i)) 



i=l 

3 



1 E( ; + o + E(- 1 )^" r E a + *i) • • • a + v) 

i=l r=2 l<ii<---<i r <j 



it follows that 

y (m - j_j = nJ y 1 

m=2 (n - m)(n - m + l)log fc m jri l(l + l)log k (n-l) 

n—1 , ., . 1 s n-2 



^(l + l)bg*(n-i) 2 ^ + log fc (n-0 



^ ^Z(/ + l)log fc (n-0' 



Applying Lemma T4. II to all the sums over I on the right hand side yields 



n-l 



£ ggj? = ^r + ££+°' 1 



m=2 



^ (n — m)(n — m + 1) log fc m Mog fc n nlog^n ^n\og k+1 n 



log fe *n log' n Mog fe+1 n'' 2 Mog fc rt Milog fc n 

A 1 rf- ] 



r=2 

i fe ife-i 
log n log n 



\og K n Mog fc+1 n 



+ (». 1 ., (T . 1) .ia+ii +D _ 1 rfO=iT)^- + o 



, _ , 7/ r- 1/ log n \log ft+1 n 

The result follows, since X^=2( _1 ) r (r) 7=1 = jhj-i -j + l = + jj - j + 1 for all j G N. □ 



The following Lemma 14.31 serves as a stopping rule for difference procedures of the form as for 
example used in the proof of Theorem 11.21 For similar stopping rules we refer the reader to [17\ 
Lemma A. 2]. 

Lemma 4.3 Fix j € N and let (a n ) ng N be a sequence of real numbers satisfying the recursion 



n-l 



a n = q n + E P nm ~7 a ™' n e {2, 3, . . .} 

m=2 ^ '3 
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for some given sequence {q n )neH- V Qn = 0(l/(nlog fe n)) for some fixed k G N, then a n — 
0(l/log k+1 n). 

Proof. By assumption there exists a constant C > such that \q n \ < C/(n log fc n) for all n G N\{1}. 
Moreover, by Corollary 14. 21 

E(m — l)j 1 n ^-4 (to — l)j 



(n)j log fc+1 m (n- l)Wj ^ (n - m)(n - m + l)\og k+1 m 

(ri + (b?+i n ~ J + °(lo e fc+1 n)) 



1 1 ^/ 1 



log fe+1 n nlog fc n V n l g fe+1 n 
Thus, there exists a constant no G N such that the inequality 



n-l 
m=2 



(m-l)j 11 i 1 



nrn 



(n)j log k+1 m - \og k+1 n 2 n \og k 



n 



holds for all integers n > no. Choose a constant Z? > 2C/j sufficiently large such that 

K\ < ? +1 (24) 
log + n 

for all n G {2, 3, . . . , no}. In the following it is verified by induction on n that ([24| holds for all 
n G N \ {1}. For n G {2, 3, ... , no}, (|2"4"|) holds by the choice of D. Suppose now that n > no and 
that \a m \ < D / log fc+1 m is already verified for all rn G {2, 3, . . . ,n — 1}. Then, 

E l (m — l)i , C ^-4 (m, — 1), Z) 



,„_, (n)j nlog fc n ^ 2 (n)j log fc+1 to 



nlog fc n Mog fc+1 n 2 n \ g k n J log fe+1 n nlog fc n \og k+1 n 



since C — jD/2 < 0. Thus, the induction is completed and the lemma is established. □ 
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